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Abstract:  The  chemisorption  of  ethylene  on  three  metal  surfaces,  Ni(l  1 1),  Pd(l  1 1)  and 
Pt(lll),  has  been  studied  by  an  approximate  molecular  orbital  method:  extended  Hiickel 
calculations  within  a  tight-binding  formalism.  The  overlap  population  between  individual 
fragment  molecular  orbital  (FMO)  of  the  adsorbate  and  each  atomic  orbital  (AO)  or  atom  of 
the  surface  can  be  projected  out  accurately.  The  relative  importance  of  each  FMO  to  the 
ethylene-metal  bonding  can  then  be  assessed.  Moreover,  interaction  diagrams  can  be 
constructed  to  describe  the  adsorbate- surface  bonding.  The  preferred  adsorption  sites 
determined  by  our  calculations  are  supported  by  the  vibrational  data  obtained  by  electron 
energy  loss  scattering.  The  activation  energy  for  the  transformation  of  a  7t- bonded  ethylene 
to  a  di-o-bonded  adsorbate  on  Pt(l  1 1)  has  been  estimated. 


Introduction 


Heterogeneous  catalysis  on  metal  surfaces  implies  enhanced  reactivity.  Facilitation 
of  a  reaction  may  occur  in  any  of  the  many  steps  of  a  surface  mechanism  —  chemisorption, 
reaction,  migration  on  the  surface,  further  reaction,  desorption.  If  we  are  to  understand 
these  reactions,  we  have  to  understand  them  in  detail  for  one  metal,  and,  most  importantly, 
how  they  vary  as  we  pass  from  one  metal  to  another. 

The  common  crystal  forms  of  the  group  10  elements  Ni,  Pd  and  Pt  are  face- 
centered  cubic.  The  easy-to-come  by  hexagonal  (1 1 1)  face  is  a  close-packed  plane,  likely 
to  be  minimally  reconstructed.  For  a  given  adsorbate,  we  might  expect  a  smooth,  gradual 

change  in  adsorbate-surface  interaction  as  we  go  from  Ni(lll)  — ►  Pd(lll)  — *  Pt( 111). 
Consider  now  specifically  ethylene.  Electron  energy  loss  scattering  (EELS)  showed  that 
ethylene  is  more  strongly  bound  and  more  distorted  on  Ni(l  1 1)  and  Pt( 111)  than  Pd(l  1 1).1 
Comparison  of  the  vibrational  spectra  also  suggested  a  different  adsorption  site  on 
Pd(l  1 1).2  This  is  inconsistent  with  the  simple  idea  of  a  gradual  change  and  poses  a  puzzle. 

In  this  paper,  we  investigate  the  chemisorption  of  ethylene  on  the  above-mentioned 
three  metal  surfaces.  The  tight-binding  extended  Hiickel  method  is  employed  to  trace  down 
the  orbital  interactions.3  As  a  compromise  between  accuracy  and  economy,  three-layer 
slabs  serve  as  models  for  the  metal  surfaces.  The  adsorbate  molecules  are  put  on  one  side 
of  the  two-dimensional  metal  slabs  only.  This  covered  layer  will  be  referred  to  hereafter  as 
the  surface  layer. 

Let  us  first  review  the  experimental  data  in  some  detail,  for  these  have  been 
subjected  to  some  ambiguity  and  controversy.  Ethylene  adsorbs  initially  in  a  disordered 
way  on  Pt(l  1 1).  After  exposure  to  the  electron  beam  (and  probably  some  decomposition)  a 
(2x2)  low  energy  electron  diffraction  (LEED)  pattern  can  be  observed,  indicating  a  possible 

coverage,  0,  of  0.25.4  This  value  of  9  was  supported  by  thermal  desorption 
spectroscopy5  and  radiotracer  data,6  but  disagreed  with  results  from  X-ray  photoelectron 


1 


spectroscopy7,8  and  elastic  recoil  detection.9  The  latter  two  methods  favored  0  =  0.5.  Our 
calculations  will  show  that  for  one-half  coverage  there  are  extremely  large  adsorbate- 
adsorbate  repulsions,  unless  the  ethylene  molecules  are  tilted  with  respect  to  the  surface. 
This  upright  position  has  been  detected  only  on  the  Pt(l  1 1)  surface  near  300  K,  at  which 
temperature  decomposition  occurred.10  We  wish  to  concentrate  here  on  low  temperature 
adsorption  where  the  C-C  bond  is  parallel  to  the  surface.  So  one-quarter  coverage  seems  to 
be  more  reasonable.  A  (2x2)  LEED  pattern  has  been  obtained  for  ethylene  on  Ni(l  1 1)  as 
well.11  The  Pt-Pt  bond  length  (2.77  A)  is  only  slightly  longer  than  that  of  Pd-Pd  (2.75  A). 
We  thus  assume  the  same  coverage,  0  =  1/4,  for  all  three  metal  surfaces. 

Although  the  exact  geometry  and  location  of  the  chemisorbed  ethylene  are 
unknown,  some  useful  information  is  available.  EELS  performed  near  100  K  pointed  to  a 

2-fold  bridge  site  (di-o-bonded  ethylene)  1  for  Ni(l  1 1)  and  Pt(  111)  but  an  on-top  site  (7t- 
bonded  ethylene)  2  for  Pd(l  1 1).2-12'15  The  assignment  of  adsorption  site  for  Ni(l  11)  was 
supported  by  near  edge  extended  X-ray  absorption  fine  structure  (NEXAFS).16  This 
technique,  however,  favored  a  rc-bonded  adsorbate  on  Pt(lll).10  Adding  to  this 
controversy  are  the  recent  results  of  ultra-violet  photoclectron  spectroscopy  (UPS)  which 
suggested  the  existence  of  a  7t-bonded  species  at  temperatures  lower  than  52  K  on  Pt(  1 1 1 ). 
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At  a  temperature  higher  than  52  K,  the  adsorbate  transforms  into  a  di-o-bonded  ethylene.17 
Finally,  it  is  interesting  to  note  that  the  di-c-bonded  ethylene  was  proposed  to  be  the 
intermediate  for  catalytic  exchange  between  etnane  and  deuterium  (the  ct(3  process)  over 
evaporated  film  of  Ni,  Pd  and  Pt  in  the  fifties  and  sixties.  More  recendy  the  involvement  of 
a  rc-bonded  alkene  has  been  put  forward.18,19  By  adjusting  the  geometry  of  the  ethylene  to 
fit  the  photoemission  results,  Demuth  suggested  a  C-C  bond  length  in  ethylene  of  1.39, 
1.44  and  1.49  A  on  Ni(lll),  Pd(lll)  and  Pt(lll)  respectively,20  (compared  to  1.37  A  in 
Zeise's  salt)  and  some  bending  of  hydrogen  away  from  the  surfaces.  The  C-C  bond  length 

on  Pt(lll)  is  in  excellent  agreement  with  the  NEXAFS  data  (1.49  ±  0.03  A).21  Room 
temperature  adsorption  of  ethylene  on  Ni(l  1 1)  produced  a  C2  surface  species  with  a  C-C 
bond  about  1.45  A.16  Recent  studies  on  Cu(100)  led  to  a  C-C  bond  lengthening  of  about 
0.13  A.  22,23 


To  simplify  the  problem  within  our  computational  ability,  all  the  C-C  bonds,  the 
C-H  bonds  and  the  HCH  angles  of  ethylene  on  Ni(lll),  Pd(lll)  and  Pt(l  11)  are  set  to 
1.45  A,  1.10  A  and  120°  respectively.  The  Ni-Ni,  Pd-Pd  and  Pt-Pt  bond  lengths  are  kept 
fixed  at  2.49,  2.75  and  2.77  A.  No  surface  reconstruction  or  relaxation  is  presumed.  From 
the  known  structures  of  organometallic  compounds  and  the  theoretical  work  of 
Anderson,24*27  we  assume  a  metal-carbon  bond  length  of  2.02  A  for  ethylene  on  Ni(l  1 1) 
and  2. 10  A  on  Pd(l  1 1)  and  Pt(l  1 1). 

The  degree  of  bending  of  the  hydrogens  should  be  different  for  each  metal  surface. 
The  bending  angles  of  the  HCH  planes  with  respect  to  the  metal  surfaces  (<j),  as  defined  in 
3),  were  determined  by  minimizing  the  total  energy.  Table  1  displays  the  optimized 
bending  angles.  Despite  our  drastic  approximations,  for  a  given  adsorption  site  ethylene 
seems  to  prefer  a  less  distorted  geometry  on  Pd(l  1 1),  in  agreement  with  the  experimental 
result. 


Table  1  here 
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Table  1.  Optimized  HCH  Bending  Angles  (<J),  see  3)  of  Ethylene  with  the  Metal  Surfaces 


3 

Before  we  proceed  to  a  detailed  analysis  of  the  adsorbate-surface  interaction,  a  few 
words  are  in  order  about  the  two  important  conceptual  tools  we  use:  the  density  of  states 
(DOS)  and  its  partitions,  and  the  crystal  orbital  overlap  population  (COOP)  curve.  The 
DOS  curve  is  a  graph  of  the  number  of  orbitals  per  unit  volume  per  unit  energy  as  abscissa, 
versus  energy  as  ordinate.  That  DOS  may  be  partitioned  on  an  atom,  atomic  orbital,  or 
fragment  molecular  orbital  basis.  The  COOP  curve  is  a  plot  of  the  overlap-population  - 
weighted  density  of  states  versus  energy.  Integration  of  the  COOP  curve  up  to  the  Fermi 
energy  then  gives  the  total  overlap  population.3  Further  computational  details  can  be 
conveniently  found  in  the  Appendix. 
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Comparison  of  the  adsorption  of  ethylene  on  Ni(lll),  Pd(lll)  and  Pt(lll) 

One  of  the  drawbacks  of  the  non-self-consistent  extended  Huckel  method  is  the 
exaggeration  of  electron  transfer.  The  valence  state  ionization  potentials  (Hu)  for  the  metals 
are  thus  refined  by  charge  iteration  of  the  three-layer  slabs  with  ethylene  at  the  2-fold  site 
and  the  on-top  site  until  self-consistency  in  the  metal  charge  distribution  is  attained.  In  this 
way,  the  excessive  electron  flow  is  greatly  reduced.  Table  2  gives  the  H-'s  of  the  s,  p  and 
d  orbitals  of  Ni,  Pd  and  Pt  so  obtained.  Among  the  three  metals,  the  energy  difference 
between  the  d  orbitals  and  the  s  or  p  is  greatest  for  palladium.  As  we  shall  see  later,  this 
leads  to  some  interesting  bonding  features.  Note  that  the  two  set  of  charge  self-consistent 
parameters  are  not  that  different 

It  is  instructive  to  compare  our  computed  results  with  those  obtained  by  other 
theoretical  methods.  Relativistic  augmented-plane-wave  calculations  suggest  that  the 
palladium  d  band  is  lower  than  that  of  the  platinum.28  This  computational  technique  also 
indicates  that  compared  to  palladium,  the  platinum  sp  bands  move  toward  lower  energies 
relative  to  the  d  band.  There  is  thus  stronger  sp-d  hybridization  in  the  case  of  platinum. 


Table  2  here 


Calculations  have  been  performed  with  both  sets  of  parameters  (those  derived  from 
the  on-top  and  2-fold  site  iterations).  Selected  results  for  the  "2-fold  site  parameters"  are 
listed  in  Table  3,  those  for  the  "on-top  site  parameters"  in  Table  4.  The  binding  energy  is 
defined  as  the  difference  in  energy  between  the  adsorbate-surface  composite  system  when 
the  planar,  undistorted  ethylene  molecule  is  removed  from  the  surface  and  when  it  is 
chemisorbed  with  the  above  specified  geometry.  A  positive  binding  energy  means  an 
attractive  interaction.  It  seems  that  on  Ni(l  1 1)  and  Pt(l  1 1)  there  is  a  slight  preference  for 
the  2-fold  site  over  the  on-top  site.  The  reverse  is  true  on  Pd(l  1 1).  Comparing  the  binding 
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Table  2.  Valence  State  Ionization  Potentials  (in  eV)  Obtained  by  Charge  Iteration  with  the 
Ethylene  at  the  2-fold  or  Atop  site. 


2-fold  site 

atop  site 

metal  surface 

s 

P 

d 

s 

P 

d 

Ni(lll) 

-7.92 

-4.18 

-11.51 

-7.79 

-4.08 

-11.30 

Pd(lll) 

-7.51 

-3.86 

-12.53 

-7.50 

-3.85 

-12.51 

Pt(lll) 

-8.82 

-5.28 

-12.15 

-8.75 

-5.23 

-12.04 

energy  at  the  2-fold  site  of  Ni(l  1 1)  and  Pt( 111)  with  that  at  the  on-top  site  of  Pd( 111) 
suggests  that  chemisorption  on  the  latter  may  be  weaker.  The  adsorption  of  ethylene  on 
Pt(l  1 1)  occurs  even  at  37  K.17  The  activation  energy  for  adsorption  should  be  vanishingly 
small.  The  activation  energy  for  desorption  would  then  be  a  good  upper  bound  for  the 
binding  energy.  The  activation  energy  for  desorption  of  the  di-a-bonded  and  the  tr-bonded 
ethylene  on  Pt(lll)  are  determined  experimentally  to  be  0.39-0.74  and  0.22-0.39  eV, 
respectively.29,30  The  computed  binding  energies  for  ethylene  on  Pt(lll)  are  probably  of 
the  right  order  of  magnitude. 

The  C-C  overlap  population  (OP)  of  ethylene  and  ethane  at  their  equilibrium 
geometries  are  1.300  and  0.737,  respectively.  Hence  the  computed  overlap  populations 
point  to  a  bond  order  around  1.5.  The  C-C  OP  for  the  more  stable  adsorption  site  on 
Ni(lll),  Pt(lll)  and  Pd(lll)  is  calculated  to  be  0.902,  0.914  and  1.066,  respectively 
(using  the  2-fold  site  parameters).  The  C-C  stretching  frequencies  for  ethylene  on  Ni(l  1 1), 
Pt(lll)  and  Pd(lll)  appear  to  be  1045,  1110  and  1355  cm1,  respectively;31  the  first  two 
values  being  the  mean  values  from  the  surface  species  derived  from  CjH4  and  C2D4,  and 
the  third  from  the  C2D4  species  alone  because  the  C2H4  spectrum  (like  that  of  Zeise's  salt 
and  of  ethylene  itself)  is  complicated  by  strong  coupling  between  C-C  stretching  and  CH2 
scissors  modes.32  An  alternate  assignment  is  1200,  1230  and  1502  cm1.1,14  Zeise'ssalt 
(KPtCl3(CH2=CH2))  and  free  ethylene  in  gas  phase  values  are  1515  and  1623  cm'1, 
respectively. 

The  agreement  between  the  computed  results  and  experiment  is  reasonable  and  the 
trend  is  not  very  parameter- sensitive  either  (compare  Table  3  and  4).  Similar  qualitative 
features  are  obtained  with  either  set  of  parameters,  except  that  the  H^'s  for  the  on-top  site 
generally  give  a  slightly  higher  binding  energy.  Hereafter  our  discussion  will  be  based  on 
the  results  derived  from  the  2-fold  site  parameters. 

Table  3  and  4  here 
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Table  3.  Binding  Energies  and  C-C  Overlap  Populations  Calculated  with  the  2-fold  Site  P 


The  numbering  of  atoms  at  the  two  adsorption  sites  of  Ni(l  1 1)  is  depicted  in  4  and 
5.  Since  (2x2)  unit  cells  have  been  chosen  (see  1  and  2),  there  is  always  a  mirror 


4  5 

plane  bisecting  and  perpendicular  to  the  C-C  bond.  Hence,  in  every  adsorption  site  there 
are  two  symmetry-equivalent  C-Ni  bonds:  the  Nij-Cj  and  Ni2-C2  bond  at  the  2-fold  site; 
the  Nij-Cj  and  Nij-Cj  bond  at  the  on-top  site.  To  compare  the  bonding  between  4  and  5, 
it  is  sufficient  to  single  out  one  of  the  two  C-Ni  bonds.  Hereafter,  whenever  we  refer  to 
the  COOP  curve  or  overlap  population  of  a  C-Ni  bond,  we  imply  the  Nij-Cj  bond,  be  it  at 
the  on-top  site  or  the  2-fold  site  (unless  otherwise  specified).  For  compactness,  the  2-fold 
site  and  on-top  site  on  Ni(lll)  are  abbreviated  as  Ni(2-fold)  and  Ni(on-top)  respectively. 
The  same  conventions  apply  to  the  other  two  surfaces. 

From  the  exp<  ience  of  organometallic  chemistry,  and  the  supporting  theoretical 
framework  of  perturbation  theory,  we  know  that  the  interactions  between  the  central  metal 
atom  and  its  ligands  occur  mainly  through  the  frontier  orbitals.  The  frontier  orbitals  of  a 
bent  ethylene  are,  in  order  of  increasing  energy,  on,  a,  n,  and  k*  6,  the  first  three  of  these 
being  fully  occupied  in  neutral  ethylene,  (o^  is  C-H  sigma-bonding  and  C-C  pi- 
antibonding.  This  is  the  inspiration  for  our  nomenclature.)  We  expect  them  to  play  an 
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important  role  on  the  adsorbate-surface  interaction.  As  we  shall  see  later,  for  the  on-top 
site  on  Pd(l  1 1)  and  the  2-fold  site  on  Ni(l  1 1)  and  Pt(  111),  these  orbitals  are  responsible 
for  125%,  99%  and  98%,  respectively  of  the  carbon-metal  OP.  (The  contribution  on 
Pd(lll)  is  greater  than  100%  because  all  the  other  eight  orbitals  of  ethylene  produce  a  net 
antibonding  OP  of  -0.041.)  Thus,  the  following  analysis  will  concentrate  on  these  four 
frontier  orbitals. 

Figure  1  is  the  analogue  of  an  interaction  diagram  in  molecular  chemistry.  In  the 
left  panel  is  the  total  density  of  states  (DOS)  of  a  monolayer  of  ethylene  arranged  in  the 
same  geometry  as  the  adsorbed  layer  at  the  2-fold  site  of  Ni(l  1 1).  The  tc  and  k*  band  then 
correspond  to  the  highest  occupied  molecular  orbital  (HOMO)  and  the  lowest  unoccupied 
molecular  orbital  (LUMO)  of  a  discrete  molecule.  There  is  extensive  overlapping  of  the 
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DOS  of  <j  and  aK  near  -15  eV.  Judging  from  the  narrowness  of  the  bands,  the  adsorbate - 
adsorbate  interaction  at  one-quarter  coverage  is  small.  This  is  expected,  for  the  closest 
H—H  contact  between  neighboring  ethylene  molecules  at  this  coverage  is  2.42  A. 
Energetically,  this  steric  interaction  costs  only  about  0.049  eV. 


Figure  1  here 


The  right  panel  is  the  DOS  of  the  bare  Ni(l  1 1)  surface.  The  position  of  the  Fermi 
energy  suggests  that  most  of  the  d  bands  are  filled.  Indeed,  for  the  surface  atoms  of  bare 
Ni(lll),  Pd(lll)  and  Pt(lll),  the  computed  electronic  configurations  are  d9  47s°-40p018, 
j9.8i5O.1OpO.09  ^  j9.36so.46p0.23j  respectively.  The  electron  occupation  of  the  valence  s 
orbitals  seems  to  be  slightly  lower,  that  of  the  p  slightly  greater  than  would  have  been 
anticipated.  The  width  of  the  d  bands  is  about  4  eV,  while  the  dispersions  of  the  s  and  p 
bands  are  much  larger,  reflecting  the  much  more  contracted  nature  of  the  d  orbitals. 

In  the  middle  panel,  we  display  the  total  density  of  states  of  ethylene  (shaded)  after 
adsorption.  In  reality,  the  Fermi  energy  of  the  metal  surfaces  should  not  move  after 
adsorption.  It  does  so  in  our  calculations,  because  of  the  finite  thickness  of  the  slabs.  The 
shifts,  however,  are  very  small  ( in  this  work  always  less  than  0. 13  eV). 

The  projected  DOS  of  the  four  important  frontier  orbitals,  magnified  by  a  factor  of 
five,  is  displayed  in  Figure  2.  The  horizontal  "sticks"  display  the  positions  of  the  MO’s  in 
a  free  ,  planar,  undistorted  ethylene  molecule.  After  adsorption,  about  18%  of  the  DOS  of 
re  is  pushed  up  above  the  Fermi  energy  and  the  main  body  of  this  band  is  pushed  down  by 
approximately  0.5  eV.  From  the  integration  curve,  about  44%  of  the  n*  gets  occupied.  It 
also  develops  substantial  dispersion,  indicative  of  strong  interaction. 
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Total  DOS  of  a  monolayer  of  ethylene  arranged  in  the  same  geometry  as  the  adsorbed 
Ni(2-fold). 

Total  DOS  of  ethylene  (shaded)  at  Ni(2-fold)  at  one-quarter  coverage. 

Total  DOS  of  Ni(l  11 )  surface. 


Figure  2  here 


Let  us  examine  the  bonding  in  these  chemisorbed  ethylenes  by  looking  at  the  crystal 
orbital  overlap  population  (COOP)  curves  in  Figure  3.  In  the  left  panel  are  the  COOP 
curves  for  the  C-Ni  bond.  There  are  three  important  bonding  regions  defined  in  the  left 
panel.  Regions  I  and  II  consist  of  the  bands  near  -14.5  eV  and  the  sharp  peak  near  -13.5 
eV,  respectively.  Region  in  starts  from  the  top  of  region  II  and  goes  up  to  the  Fermi 

energy.  Over  90%  of  the  DOS  of  o  and  cr,,.  can  be  found  in  the  first  region,  while  the  it 
band  is  the  major  component  of  region  II.  A  substantial  amount  of  the  DOS  of  7t  (-10%), 
and  nearly  all  the  DOS  of  it*  that  is  below  the  Fermi  energy,  is  located  in  region  HI.  The 
bonding  in  region  n  is  almost  entirely  due  to  it-to- metal  forward  donation,  while  in  region 
HI  it  is  mainly  the  result  of  metal-to-it*  back-bonding  and  partly  that  of  forward  donation 
involving  it. 


Figure  3  here 


A  similar  distribution  of  the  DOS  of  the  four  frontier  orbitals  among  the  three 
bonding  regions  occurs  at  the  on-top  site  of  Ni(lll)  and  the  other  two  metal  surfaces, 
Pd(l  11)  and  Pt(l  1 1).  There  are,  however,  two  major  trends.  The  DOS  of  it*  in  region  III 
decreases  from  the  2-fold  site  to  the  on-top  site.  For  a  given  adsorption  site,  it  increases  in 
the  following  order:  Pd(lll)  — » Pt( 111)  — »  Ni(lll). 

We  can  partition  the  contributions  to  adsorbate-surface  bonding  by  energy  region; 
this  and  other  bonding  information  is  listed  in  Table  5.  From  the  electron  population  of  it 
and  it*,  the  strongest  forward  donation  and  back-bonding  seem  to  occur  at  the  2-fold  site 
of  Ni(l  11)  and  Pt(  1 1 1 ),  respectively. 
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(A®)  A6i»u3 


The  solid  line  and  the  dotted  line 


Figure  3.  a)  COOP  curve  for  C-Ni  bond. 

b)  COOP  curve  for  C-Pd  bond. 

c)  COOP  curve  for  C-Pt  bond. 


Table  5  here 


The  change  in  dispersion  and  electron  occupation  of  o  and  oK  is  not  dramatic. 
Does  this  mean  that  c  and  <Jn  do  not  participate  significantly  in  bonding  to  the  surface  ? 
Although  a  and  an  show  little  change  in  electron  occupation  and  do  not  develop  substantial 
dispersion,  they  do  interact  with  the  surface,  especially  at  the  on-top  site.  For  example, 
region  I  accounts  for  about  17%  of  the  C-Pd  OP  at  Pd(2-fold)  but  44%  at  Pd(on-top).  The 
COOP  curves  of  Pt(on-top)  in  Figure  3  shows  that  the  bonding  peak  in  region  I  is  actually 
larger  than  that  in  region  II. 

Our  next  goal  is  to  determine  whether  o  or  an  is  the  major  source  of  ethylene-metal 
bonding  in  this  region.  From  its  shape,  o  should  overlap  much  better  with  the  metal 
orbitals  at  the  on-top  site,  as  shown  in  7.  The  reverse  should  be  true  for  <JK,  because  at 


7 

the  on-top  site  it  is  locally  orthogonal  to  all  orbitals  of  the  nearest  metal  atoms.  From  the 
OP  in  Table  5  and  the  COOP  curves  in  Figure  3,  the  ethylene-metal  interaction  in  region  I 
decreases  on  going  from  the  on-top  to  the  2-fold  site,  suggesting  the  dominance  of  a  over 
o^.  This  was  confirmed  by  a  numerical  experiment,  in  which  selected  atomic  overlaps 
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Table  5.  Selected  Bonding  Information  at  the  2-fold  site  and  the  atop  site. 


between  the  surface  and  the  adsorbate  (those  involved  in  oK)  were  dropped. 

Since  a  is  quite  important  at  the  on- top  site  of  Pd(l  1 1)  (in  region  I,  it  accounts  for 
42%  of  the  total  C-Pd  OP),  it  is  worthwhile  to  trace  down  the  metal  orbital  that  interacts 
strongly  with  it.  Considering  its  orientation,  a  good  candidate  is  dz2  (s  and  pz  are 

eliminated  because  they  would  lead  to  substantial  depopulation  of  cr).  As  shown  in  Figure 
4,  region  I  picks  up  about  8%  of  the  DOS  of  nearest  surface  dp  orbital  at  Pd(on-top). 
Even  in  the  absence  of  the  adsorbate,  the  s,  pz  and  dp  of  the  surface  atoms  have  the  same 
symmetry  and  extensive  hybridization  is  possible.  It  is  therefore  not  a  good  idea  to 
estimate  the  OP  due  to  dp  by  dropping  the  overlap.  Instead  we  project  out  its  contribution. 
In  this  way,  mixing  of  orbitals  induces  no  error.  The  projected  OP  for  dp  is  in  region  I  is 

0.043.  The  total  projected  OP  for  all  other  metal  orbitals  is  0.030.  The  interaction  of  a 
with  the  dp  in  this  region  is  thus  significantly  greater  than  that  of  the  total  of  the  other  eight 
metal  orbitals. 


Figure  4  here 


7t  and  7t*  are  the  most  active  participants  in  the  adsorbate-surface  interaction  in 
region  II  and  HI.  They  undergo  significant  depopulation  and  population,  respectively, 
upon  adsorption.  Their  electron  occupations  also  determine  the  net  charge  on  the 
adsorbate,  which  is  roughly  equal  to  two  minus  their  total  electron  population  (see  Table 
5).  On  Ni(lll)  and  Pt(lll),  the  metal-to-7t*  back-donation  of  electron  outweighs  the 
forward  donation  from  k  and  the  adsorbed  ethylene  is  negatively  charged.  The  opposite 
occurs  on  Pd(l  11).  As  mentioned  before,  this  may  be  due  to  its  lower  d  bands.  Notice 
that  the  o  band  is  always  below  the  d  block.  A  lower  for  the  d  orbitals  means  a  better 
energy  match  for  interaction.  The  percentage  contribution  of  a  in  region  I  to  the  metal- 
carbon  OP  is  indeed  greatest  at  the  on-top  site  of  Pd(l  1 1). 

Detailed  analysis  of  region  II  and  ID,  however,  is  not  so  simple  because  it  involves 
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Proj  DOS 
Integration 


contributions  from  o,  n  and  it*,  some  bonding,  some  antibonding.  The  analytical  tools  we 
have  used  so  far  may  not  give  us  a  sufficiently  detailed  picture  of  adsorbate-surface 
interaction.  What  we  need  is  a  method  to  project  out  the  individual  contribution  from  each 
FMO  of  ethylene  to  the  carbon-metal  OP.  In  Table  6,  the  result  of  partitioning  of  the  total 
carbon-metal  overlap  populations  by  this  method  is  presented. 


Table  6  here 


Table  6  reveals  several  bonding  features.  First,  we  recall  that  there  is  strong 
interaction  between  a  and  the  metal  surfaces  at  the  on-top  site.  The  cr-metal  OP's, 
however,  are  not  significant.  In  interpreting  this  apparent  paradox,  the  o-Pd  COOP  curve 
of  Pd(on-top)  in  Figure  5  is  instructive.  The  bonding  peak  in  region  I  is  essentially 
neutralized  by  antibonding  interaction  in  region  II  and  III.  Although  there  are  good 
overlaps  between  o  and  the  d  orbitals,  the  net  effect  is  almost  nonbonding.  The  fact  that 
there  is  a  substantial  amount  of  DOS  of  <^2  in  resonance  with  that  of  o  only  means  strong 
interaction  but  not  necessarily  strong  bonding.  To  have  strong  bonding,  the  antibonding 
counterpart  has  to  be  above  the  Fermi  energy.  This  is  exactly  what  happens  when  7t  and 
it*  interact  with  the  metal.  An  analysis  of  bonding  that  is  solely  based  on  the  DOS  thus 
may  be  insufficient. 


Figure  5  here 


On  Ni (111),  probably  because  of  the  higher  Fermi  energy,  k*  always  contributes 
more  to  the  metal -carbon  overlap  population  than  7t.  The  opposite  occurs  on  Pd(l  1 1).  For 
Pt(l  1 1),  a  transition  can  be  observed:  at  the  on-top  site,  k  is  more  important,  but  at  the  2- 
fold  site  7t*  is  the  major  source  of  C-Pt  bonding.  According  to  the  electron  populations  in 
Table  5,  the  COOP  curves  in  Figure  3  and  the  projected  OP,  the  interaction  of  n  and  it* 
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Table  6.  Contribution  of  the  four  FMO's  of  Ethylene  to  the  Carbon-metal  Overlap 
Populations.  The  contributions  of  the  other  eight  FMO’s  are  negligible. 


2-fold  site 

atop  site 

FMO 

Ni(lll) 

Pd(lll) 

Pt(lll) 

Ni(lll) 

Pd(lll) 

Pt(lll) 

71* 

0.236 

0.108 

0.215 

0.118 

0.072 

0.113 

7t 

0.137 

0.119 

0.172 

0.108 

0.118 

0.140 

a 

-0.001 

-0.003 

0.001 

0.032 

0.018 

0.043 

-0.005 

-0.002 

-0.003 

-0.002 

-0.001 

-0.002 

Total  metal-carbon 

overlap  population 

0.371 

0.212 

0.394 

0.229 

0.165 

0.271 

-7 


-9H 


-iH 


Figure  5.  a- Pd!  COOP  curve. 


with  the  metal  surface  strengthens  on  going  from  the  on-top  site  to  the  2-fold  site.  This 
may  be  a  combination  of  two  factors:  orientation  and  symmetry.  The  lobes  of  7t  and  K*  are 
pointing  towards  the  metal  atoms  at  the  2-fold  site  but  away  from  them  at  the  on-top  site,  8. 


8 

For  a  Jt-bonded  ethylene,  symmetry  prevents  any  good  overlap  between  n  and  the 
nearest  metal  px  and  dxz.  It  also  prohibits  any  effective  interaction  of  7t*  with  the  surface  s, 
pz  and  dz2  9.  There  is  another  interesting  but  more  subtle  feature  concerning  rc*.  For 
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Ni(lll)  and  Pt(lll),  the  and  OP  decrease  significantly  (  50%  for  the 

former  and  46%  for  the  latter),  on  changing  :'rom  a  di-ct-bonded  ethylene  to  a  7t-bonded 
adsorbate.  The  corresponding  reduction  of  OP  on  Pd(lll)  is  much  smaller,  only  about 
8%.  In  addition,  if  we  compare  the  electron  occupation  of  it*  between  the  two  adsorption 
sites  among  the  three  metal  surfaces,  again  there  is  a  much  smaller  difference  on  Pd(l  1 1). 

This  differential  weakening  of  7t*-metai  interaction  may  be  understood  by 
examining  Figure  6.  In  the  left  panel  are  the  COOP  curves  for  7t*-Pd1  at  Pd(2-fold)  and 
Pd(on-top).  By  projecting  out  the  OP  between  it*  and  individual  palladium  orbitals,  it  is 
found  that  the  two  bonding  peaks  at  -12.5  eV  (within  region  IQ)  are  mainly  the  result  of 
interaction  between  it*  and  d^  at  the  on-top  site;  it*,  and  dxz  at  the  2-fold  site.  These 
two  peaks  are  of  similar  size  and  hence  the  similar  7t*-Pdj  OP  at  Pd(2-fold)  and  Pd(on- 
top).  Near  -9.5  eV,  where  most  of  the  DOS  of  it*  is  located,  the  situation  is  completely 
different.  For  Pd(on-top),  there  is  some  rc*-px  bonding  interaction  which  is  almost  exactly 
cancelled  out  by  the  antibonding  interaction  between  n*  and  d^.  The  net  result  is  that  the 
orbitals  in  this  region  are  mainly  7t*-Pd  nonbonding.  For  Pt(2-fold),  the  bands  near  -9.5 
eV  are  ?t*-dxz  and  it*-^  antibonding  but  Jt*-s,  Jt*-pz  and  it*-px  bonding.  The  bonding 
interactions  of  s  and  p  outweigh  the  antibonding  interactions  of  the  d  and  a  prominent 
bonding  peak  appears.  ( This  peak  actually  is  larger  than  that  in  region  m). 


Figure  6  here 


In  summary,  the  7t*-Pd  bonding  at  the  on-top  site  may  be  represented  schematically 
by  a  'three-level'  interaction  diagram  in  10.  A  characteristic  outcome  of  this  type  of 
interaction  is  a  nonbonding  level  in  the  middle.  In  this  case,  it  corresponds  to  the  it* 
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(A*)  *6i»u3 


-Ni  COOP  curve. 


orbitals  just  above  the  Fermi  energy.  At  the  2-fold  site,  if  we  consider  px  and  py  as  one 
level,  dxz  and  dz2  as  another,  a  'four-level'  interaction  diagram  of  the  type  shown  in  11 
may  be  appropriate  to  describe  the  rt*-Pd  bonding.  The  second  level  ,which  is  labelled  as 
A,  is  the  cause  of  the  large  bonding  peak  near  -9.5  eV.  The  orbitals  on  the  left  side  of  this 
drawing  are  meant  to  remind  the  reader  that  re*  can  overlap  effectively  only  with  those 
metal  orbitals  that  are  Pdj-Pt^  antibonding. 


The  JC*-Ptj  COOP  curves  are  displayed  in  the  middle  panel  of  Figure  6.  Part  of  A 
has  been  pushed  down  below  the  Fermi  energy  and  the  bonding  peak  near  -9.5  eV  is 
smaller  than  that  in  region  HI.  The  result  is  dramatic.  In  Pd(2-fold),  the  s  and  p  bands 
only  account  for  27%  of  the  rt*-Pd  OP  but  in  Pt(2-foId),  their  contribution  to  the  7t*-Pt  OP 
is  about  50%.  The  K^-surface  bonding  in  the  upper  part  of  region  HI  is  dominated  by  the 
s  and  p  bands.  For  Ni  (1 1 1),  as  indicated  by  the  COOP  curves  in  the  right  panel,  an  even 
larger  portion  of  A  has  penetrated  into  region  III.  Now  over  59%  of  the  Jt*-Ni  OP  can  be 
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attributed  to  the  s  and  the  p  orbitals.  These  results  also  fit  the  general  trend  that  the 
chemistry  of  palladium  is  more  similar  to  that  of  platinum  than  nickel.26 

Thus,  it  is  the  position  of  A  relative  to  the  Fermi  energy  that  is  mainly  responsible 
for  the  differential  reduction  in  the  OP  between  n*  and  the  nearest  surface  atom  on  going 
from  the  2-fold  to  the  on-top  site.  The  incapability  of  the  s  and  p  orbitals  of  the  palladium 
to  push  A  down  below  the  Fermi  energy  greatly  diminishes  the  back-bonding,  and  Pd(2- 
fold)  is  the  only  2-fold  site  at  which  n  contributes  more  to  the  carbon-metal  OP  than  7t*. 

A  more  thorough,  quantitative  description  of  the  bonding  of  K  and  n*  with  the 
metal  orbitals  is  given  in  Table  7  and  8.  Due  to  a  symmetry  constraint,  almost  half  of  the 
OP  in  these  two  tables  are  equal  or  close  to  zero.  The  OP's  in  Table  7  clearly  show  that  the 
re-surface  bonding  occurs  mainly  via  the  metal  s  and  p  bands,  in  agreement  with  the 
concept  of  forward  donation.  Similarly,  the  idea  of  back-bonding  is  well  illustrated  by  the 
results  in  Table  8. 


Table  7  and  8  here 


So  far,  we  have  concentrated  on  the  adsorbate.  How  about  the  surfaces  ? 
Generally,  the  greater  the  number  of  surface  atoms  directly  bonded  to  the  adsorbate,  the 
greater  will  be  the  bond- weakening  within  the  surface.  The  metal-metal  OP's,  however, 
indicate  that  the  degree  of  bond  weakening  at  the  surface  for  the  on-top  and  2-fold  site  is 
essentially  the  same,  that  is,  this  factor  is  not  crucial  in  determining  the  site  preference.  As 
we  shall  see  in  a  future  paper,  such  is  not  the  case  if  the  adsorbate  is  carbon  monoxide. 
One  would  like  to  separate  the  carbon-metal  OP  into  nine  components,  corresponding  to  the 
contribution  of  each  orbital  of  the  metal  atoms  involved  in  anchoring  the  ethylene.  The  s 
and  p  bands  of  a  transition  metal  are  very  diffuse.  Their  band-widths  are  tens  of  eV.  The 
changes  in  DOS  of  these  orbitals  after  adsorption  are  thus  much  less  obvious  than  those  of 
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Table  7.  Overlap  Population  between  7t  and  the  Nine  Metal  Surface  Orbitals.  The  py,  dxy, 
d^  interactions  vanish  by  local  symmetry,  and  are  not  listed. 


2-fold  site 

atop  site 

metal  orbital 

Ni(lll) 

Pd(lll) 

Pt(lll) 

Ni(lll) 

Pd(lll) 

Pt(lll) 

s 

0.075 

0.082 

0.082 

0.067 

0.085 

0.077 

Px 

0.001 

0.002 

0.003 

0 

0 

0 

Pz 

0.042 

0.031 

0.050 

0.0290 

0.027 

0.039 

dxV 

0 

0 

0 

-0.001 

-0.001 

-0.001 

dz2 

0.016 

0.004 

0.028 

0.013 

0.006 

0.025 

dxz 

0.003 

0.001 

0.009 

0 

0 

0 

Table  8.  Overlap  Population  between  k*  and  the  Nine  Metal  Surface  Orbitals.  The  py, 
dxy,  dyZ  interactions  vanish  by  local  symmetry,  and  arc  not  listed. 


the  d  bands.  As  a  result,  the  role  of  the  s  and  p  orbitals  on  chemisorption  involving 
transition  metals  is  seldom  discussed,  although,  by  analogy  to  organometallic  compounds, 
these  orbitals  should  interact  with  the  adsorbate.  Their  counterparts,  the  more  contracted  d 
orbitals,  may  be  over-emphasized  in  the  existing  literature  discussions,  including  our  own. 
The  partitioning  of  the  carbon-metal  OP  (Table  9)  may  provide  more  insight  on  this  issue. 


Table  9  here 


There  is  little  doubt  that  the  d  orbitals  are  involved  in  adsorbate-surface  bonding, 
but  the  combined  effect  of  the  s  and  p  bands  is  much  more  important.  In  fact,  apart  from 
the  on-top  site  on  Ni(l  1 1),  the  s  orbital  alone  contributes  more  to  the  carbon-metal  OP  than 
the  total  of  the  five  d  orbitals.  As  shown  by  the  OP  in  Table  7  and  8,  the  reason  behind  this 
is  that,  apart  from  forward  donation,  the  surface  s  and  p  bands  participate  actively  in  k*- 

metal  bonding.  In  all  adsorption  sites,  py,  dxy  and  dyz  are  orthogonal  to  a,  7t  and  n*. 
Most  of  the  interaction  between  px  and  metal  surfaces  occurs  above  the  Fermi  energy. 
&,2.y2  may  engage  heavily  in  metal-metal  bonding.  Hence  these  five  orbitals  only  play  a 
minor  role  in  surface-adsorbate  bonding. 

The  interaction  between  the  ethylene  and  the  surface  s,  pz  and  6^2  bands  reduces 
substantially  as  one  moves  from  the  2-fold  to  the  on-top  site.  From  the  previous 
discussion,  we  think  that  this  may  be  related  to  the  symmetry  of  ji*.  To  confirm  our 
conjecture,  we  take  a  close  look  at  the  COOP  curves  between  Ni  s  and  carbon  in  Figure  7. 
For  the  2-fold  site,  there  is  significant  bonding  contribution  from  region  III,  indicative  of 
strong  interaction  with  the  7t*.  k*  is  orthogonal  to  the  nearest  surface  s  orbital  at  the  on-top 
site.  Effective  overlap  is  thus  prohibited  and  the  bonding  interaction  is  much  smaller.  A 
similar  argument  seems  to  apply  to  the  pz  and  dj,2. 
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Table  9.  Contribution  of  the  Nine  Metal  Orbitals  to  the  Carbon-metal  Overlap 
Populations. 


2-fold  site _ _ atop  site 


metal  orbital 

Ni(lll) 

Pd(lll) 

Pt(lll) 

Ni(lll) 

Pd(lll) 

Pt(lll) 

s 

0.169 

0.117 

0.160 

0.095 

0.105 

0.105 

Px 

0.004 

0.005 

0.010 

0.019 

0.007 

0.019 

Py 

-0.001 

-0.000 

0.001 

-0.002 

-0.003 

-0.001 

Pz 

0.092 

0.015 

0.089 

0.022 

-0.000 

0.036 

dxV 

-0.002 

-0.002 

-0.002 

-0.001 

-0.001 

-0.001 

V 

0.104 

0.062 

0.110 

-0.003 

-0.001 

0.018 

dxy 

0.001 

-0.000 

0.000 

-0.001 

-0.0010 

-0.001 

dxz 

0.010 

0.020 

0.031 

0.101 

0.067 

0.097 

dyz 

-0.005 

-0.005 

-0.005 

-0.000 

-0.001 

-0.001 

Figure  7  here 


In  contrast  to  other  metal  orbitals,  the  projected  OP  of  is  much  larger  at  the  on- 
top  site.  To  trace  down  the  reason  for  this,  we  examine  the  COOP  curves  for  the  carbon 
atom  and  the  Ni  d^  in  the  left  panel  of  Figure  8.  At  the  on-top  site,  d^  is  orthogonal  to  a, 
an  and  Jt  of  ethylene.  The  interactions  in  the  first  two  regions  are  thus  negligible.  For  the 
2-fold  site,  the  symmetry  constraint  is  partially  relaxed,  o,  K  and  the  pan  of  the  surface  d^ 
band  that  is  Ni}-Ni2  bonding  undergo  a  ’three-level'  interaction.  However,  the  greatest 
difference  occurs  in  region  m.  In  Ni(2-fold),  it*  can  only  overlap  effectively  with  the 
upper  part  of  the  surface  d^  band  that  is  Nij-Ni2  antibonding.  In  Ni(on-top),  it  can 
interact  extensively  with  the  whole  d^  band,  as  shown  in  12.  Thus  a  much  larger  C-dX2 


12 

OP  results  (see  Table  8).  The  COOP  curve  in  the  right  panel  is  meant  to  illustrate  the 
characteristics  of  the  above  'three-lever  interaction;  the  lowest  level  is  strongly  bonding, 
the  middle  level  essentially  nonbonding,  and  the  third  level  strongly  antibonding. 
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Figure  7.  Contribution  of  the  nickel  s  orbital  to  the  Ni-C  COOP  curve. 


Figure  8  here 


It  is  instructive  to  examine  the  partitioning  of  the  Ni-C  COOP  curves  into  the  s,  p 
and  d  components.  As  expected,  in  Figure  9  the  s  and  p  bands  produce  large  bonding 
peaks  in  region  II,  where  most  DOS  of  7t  can  be  found,  corresponding  the  rc-to-metal 
forward  donation.  Surprisingly,  this  is  also  true  for  the  d  orbitals.  This  would  induce 
some  antibonding  interaction  within  region  HI.  Hence  the  unexpectedly  small  contribution 
of  the  d  bands  to  the  carbon-metal  OP  (relative  to  s  and  p)  is  due  to  the  antibonding  peaks 
in  region  HI,  as  a  direct  result  of  interactions  of  the  d  bands  with  filled  orbitals  of  ethylene. 
(We  see  here  a  surface  analogue  of  the  molecular  two-level-four  electron  repulsion).3  The 
partitioning  of  the  C-Pd  and  C-Pt  COOP  curves  shows  similar  features. 

Figure  9  here 
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Figure  8.  a)  COOP  curve  for  the  dxz  orbital  of  surface  nickel  atoms  and  ethylene. 

b)  COOP  curve  for  the  interaction  between  the  0  and  n  orbital  of  ethylene 
with  the  d„..  orbitals  of  the  surface  nickel  atoms. 


The  Organometallic  Analogy 


It  is  instructive  to  compare  the  computational  results  for  the  three  on-top  sites  with 

those  of  discrete,  molecular  T}2  coordinated  (^-bonded)  olefinic  complexes,  studied  in  a 
previous  paper  of  our  group.33  Unfortunately  the  ethylenes  were  kept  planar  in  the 
molecular  calculations,  for  reasons  of  simplicity.  For  comparison,  we  planarize  our 
adsorbates  and  redo  the  calculations.  Selected  information  is  collected  in  Table  10,  where 
the  occupation  of  n  and  K*  of  ethylene  for  molecules  13  -  18  is  reported.  Note  that  the 


generally  similar  occupations  for  discrete  molecules  and  surfaces  for  a  given  metal,  and  the 
wider  range  of  occupations  in  known  organometallic  olefin  complexes. 


Table  10  here 
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Table  10.  Electron  Density  of  n  and  K*  and  the  Net  Charge  on  Ethylene. 


System 

Jt 

n* 

Net  Charge 

Ni(PH3)2C2H4 

1.81 

0.44 

-0.28 

Ni(on-top) 

1.76 

0.47 

-0.27 

(PdCl3C2H4)- 

1.71 

0.21 

+0.02 

Pd(on-top) 

1.66 

0.35 

+0.01 

cis-PtCl2NH3C2H4 

1.64 

0.24 

+0.10 

Pt(on-top) 

1.66 

0.35 

+0.01 

CpFe(CO)2C2H4+ 

1.53 

0.32 

+0.20 

Fe(CO)4C2H4 

1.66 

0.48 

-0.15 

Cp2WHC2H4+ 

1.74 

0.60 

-0.36 

There  are  certainly  many  similarities  between  the  organometallic  compounds  and 
surface  complexes.  Yet  there  are  differences  too.  For  instance,  the  energy  of  the  HOMO 
depends  strongly  on  the  nature  of  the  ligand,  while  the  metal  Fermi  energy  is  not  affected 
by  the  adsorbate  molecules.  The  work  function  of  the  metal  does  vary  after  chemisorption, 
due  to  a  change  in  the  surface  dipole.  The  similarities  build  the  connections.  The 
differences  make  each  of  them  an  exciting  research  subject  on  its  own. 
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Interconversion  of  ^-bonded  to  di-a-bonded  ethylene 

To  our  knowledge,  there  have  been  two  reports  on  the  transformation  of  a  71- 

bonded  ethylene  to  a  di-a-bonded  adsorbate  on  Pt(l  11). 17,34  We  would  like  to  estimate 
the  activation  energy  of  this  process.  Given  our  geometries,  the  transformation  involves  a 

translational  motion  of  1.3852  A  and  an  increase  in  the  bending  angle  of  the  HCH  plane  (<j>, 
see  3)  by  14°.  We  thus  assume  a  reaction  pathway  in  which  there  is  an  increase  in  1.4°  of 
bending  angle  for  every  linear  translation  of  0.13852  A,  19.  In  Figure  10,  we  present  the 
potential  energy  curve  for  this  transformation.  The  activation  energy  was  found  to  be 
0.22  eV.  Hence  the  transformation  seem  to  be  energetically  feasible,  this  is  a  small  energy. 


<f>  =  39° 
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By  assuming  a  similar  reaction  pathway  (one  degree  decrease  in  bending  angle  0  per  linear 
translation  of  1.3752  A),  the  activation  energy  for  the  conversion  of  the  2-fold  site  on 
Pd(  111)  to  the  on-top  site  was  estimated  to  be  0.14  eV. 
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There  have  been  some  indications  of  the  existence  of  a  rc-bonded  ethylene  on 
Pt(l  1 1),12,17,34  Since  the  energy  differences  for  the  two  adsorption  sites  on  the  two 
surfaces,  Pd(l  11)  and  Pt(l  1 1)  are  similar  (but  in  opposite  order)  and  the  activation  energies 

for  transformation  are  roughly  the  same,  it  may  be  worthwhile  to  try  to  detect  the  di-o- 
bonded  species  on  Pd(l  1 1).  On  the  contrary,  the  energy  barrier  for  the  transformation  of 
ethylene  from  the  on-top  site  on  Ni(lll)  to  the  2-fold  site  was  calculated  to  be  a  tiny 
0.0026  eV.  The  rc-bonded  species  on  Ni(lll)  seems  to  be  thermodynamically  and 

kinetically  unstable  with  respect  to  the  di-c-bonded  ethylene  and  its  detection,  may  be  more 
difficult 


Figure  10  here 
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Figure  10.  Energy  profile  for  the  interconversion  of  ethylene  between  Pt(2-fold)  and  Ft(on-top). 


Conclusion 


We  have  investigated  in  detail  the  adsorption  of  ethylene  on  three  metal  surfaces, 
Ni(lll),  Pd(lll)  and  Pt(lll).  Several  established  and  useful  concepts  in  molecular 
chemistry,  such  as  the  Dewar-Chatt-Duncanson  model,  local  symmetry,  three-level 
interactions  are  valuable  for  understanding  the  adsorbate-surface  interaction.  There  are, 
however,  some  unexpected  results  as  well.  It  seems  that  besides  the  d  bands,  the  s  and  p 
bands  also  play  a  significant  role  in  the  bonding  between  n*  and  the  surface.  Although  the 
o  orbital  of  ethylene  interacts  strongly  with  the  metal  atoms,  the  contribution  of  this  orbital 
to  bonding  is  not  significant  because  most  of  the  a-surface  antibonding  states  are  below  the 
Fermi  energy.  In  contrast,  most  of  the  7t-surface  and  Jt*-surface  bonding  states  are 
occupied  while  their  antibonding  counterparts  are  empty.  Hence,  these  two  frontier  orbitals 
dominate  the  adsorbate-surface  bonding.  n*  is  the  major  contributor  to  the  carbon-metal 
OP  for  Ni(2-fold)  ,Ni(on-top),  and  Pt(2-fold),  but  n  is  more  important  for  the  other  three 
adsorption  sites.  On  Ni(  1 1 1 )  and  Pt(  111),  the  Jt*-metal  bonding  is  much  stronger  at  the  2- 
fold  site  than  the  on-top  site,  mainly  due  to  level  A  (  see  scheme  11).  This  may  be  related 
to  the  slight  energetic  preference  for  Ni(2-fold)  and  Pt(2-fold)  relative  to  Ni(on-top)  and 
Pt(on-top),  although  at  the  2-fold  site  there  is  more  weakening  of  the  C-C  bond.  For  Pd(2- 
fold),  most  of  the  rc*-s  and  7t*-p  bonding  states  are  above  the  Fermi  energy.  As  reflected 
by  the  electron  occupation  and  the  projected  OP,  the  7t*-Pd  interaction  is  only  slightly 
strengthened  on  going  from  Pd(on-top)  to  Pd(2-fold).  This  small  gain  in  bonding  may  not 
be  enough  to  compensate  for  the  C-C  weakening,  and  thus  Pd(on-top)  is  the  more  stable 
site. 


26 


Appendix 


Our  calculations  are  of  the  tight-binding  extended  Hiickel  type  with  a  weighted  HtJ 
approximations.35,36  The  geometrical  assumptions  concerning  bond  lengths  (in  A) 
include  the  following:  C-H  1.10,  C-C  1.45,  Ni-C  2.02,  Pd-C  2.10,  Pt-C  2.10,  Ni-Ni 
2.49,  Pd-Pd  2.75  and  Pt-Pt  2.77  A.  The  extended  Hiickel  parameters  are  listed  in  Table 
10.  Double  zeta  expansions  of  the  metal  d  orbitals  have  been  employed.  The  parameters  of 
ethylene  are  taken  from  previous  study.37  Charge  iterations  have  been  performed, 
assuming  a  quadratic  dependence  of  metal  Hj^'s  on  charge.38 

The  detailed  procedure  for  obtaining  the  metal  parameters  is  as  follows:  First  we 
assume  a  bending  angle  of  30°  for  all  six  adsorption  sites.  We  performed  the  charge 
iterations  to  get  a  set  of  H^'s,  calling  them  A.  We  then  optimized  the  bending  angles  0, 
with  set  A.  Redoing  the  charge  iterations  with  these  0,  we  obtained  another  set  of 
parameters,  set  B.  Since  set  A  was  almost  equal  to  set  B,  0  remained  unchanged  even  if 
we  repeated  the  optimization  of  bending  angles  with  set  B.  The  parameters  used  and  listed 
in  this  paper  are  set  B. 


Table  1 1  here 


The  k  points  are  generated  according  to  the  geometrical  method  of  Btihm  and 
Ramirez.39  Since  all  the  six  adsorption  sites  have  the  same  symmetry,  the  same  k  point  set 
(totally  18)  can  be  and  has  been  chosen  for  all  calculations.  In  this  way,  the  error  in 
computing  the  energy  difference  between  the  on-top  and  2-fold  site  should  be  minimized. 
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Table  1 1 .  Extended  Hiickel  Parameters. 


Atom 

Orbital 

Hu(eV) 

Ci 

C2 

cl 

c2 

C 

2s 

-21.4 

1.62 

2p 

-11.4 

1.62 

H 

Is 

-13.6 

1.30 

Ni 

4s 

-7.92 

2.10 

4p 

l 

OO 

2.10 

3d 

-11.51 

5.75 

2.00 

0.57 

0.63 

Pd 

5s 

-7.51 

2.19 

5p 

-3.86 

2.15 

4d 

-12.53 

5.98 

2.61 

0.55 

0.67 

Pt 

6s 

-8.82 

2.55 

6p 

-5.28 

2.55 

5d 

-12.15 

6.01 

2.70 

0.63 

0.55 
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